Sums over geometries and improvements on the mean field approximation 
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The saddle points of a Lagrangian due to Efetov are analyzed. This Lagrangian was originally 
proposed as a tool for calculating systematic corrections to the Bethe approximation, a mean-field 
approximation which is important in statistical mechanics, glasses, coding theory, and combinatorial 
optimization. Detailed analysis shows that the trivial saddle point generates a sum over geometries 
reminiscent of dynamically triangulated quantum gravity, which suggests new possibilities to design 
sums over geometries for the specific purpose of obtaining improved mean field approximations to D- 
dimensional theories. In the case of the Efetov theory, the dominant geometries are locally tree-like, 
and the sum over geometries diverges in a way that is similar to quantum gravity's divergence when 
all topologies are included. Expertise from the field of dynamically triangulated quantum gravity 
about sums over geometries may be able to remedy these defects and fulfill the Efetov theory's 
original promise. The other saddle points of the Efetov Lagrangian are also analyzed; the Hessian 
at these points is nonnormal and pseudo-Hermitian, which is unusual for bosonic theories. The 
standard formula for Gaussian integrals is generalized to nonnormal kernels. 
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I. INTRODUCTION 



Modern physics often struggles with geometry. Zero 
dimensional and "infinite dimensional" systems are the 
simplest to analyze, precisely because in such systems 
the effects of dimension and distance are grossly simpli- 
fied. In this paper we will restrict our attention to the 
more realistic and challenging scenario: a finite number 
of dimensions D. In this scenario physicists typically 
choose ground states - saddle points - that are spatially 
uniform, which means that we are pretending that the 
system "almost" lacks spatial structure. Geometrical in- 
formation is later added back into perturbative calcu- 
lations via the bare Green's function. This program of 
mean field theory fails when the number of dimensions is 
small, and in disordered and complex systems often must 



be augmented with replicas and replica symmetry break- 
ing [l| . It would be desirable to find a systematic pertur- 
bative expansion which incorporates geometrical effects 
non-perturbatively; in the zeroth-order approximation. 
Such a theory would include incomplete geometrical in- 
formation in the lowest order, and would systematically 
improve its geometrical accuracy at higher orders. 

One very tempting possibility is to design an expan- 
sion whose zeroth order is the Bethe approximation [24j . 
Bethe required that if a spin residing at node a is re- 
moved then the spins on neighbors of a will be uncorre- 
cted @,H|. On a tree his approximation is exact; thus its 
physical meaning consists of including complete informa- 
tion about nearest neighbors, but leaving out the effects 
of loops, i.e. possibilities to make round trips without 
retracing one's steps. When treating low-dimensional, 
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disordered, or frustrated systems, the Bethe approxima- 
tion gives notable improvements over mean field theory 
It is equivalent to using the replica method without 
breaking replica symmetry, and has a hierarchy of simple 
generalizations which are equivalent to successive levels 
of replica symmetry breaking 

A systematic expansion of corrections to the Bethe 
approximation would adjust for the difference between 
the correct D-dimensional geometry and a tree. Several 
authors have proposed improvement schemes based on 
higher correlated moments [4[ , effects of nearby spins on 
one another [j| @ , or ''generalized loops" moving around 
the correct geometry 0, HI ■ All of these schemes involve 
an exponentially large number of corrections, as is typi- 
cal in perturbative calculations. Less typically, in no for- 
malism has anyone identified a small parameter whose 
powers can be used to justify truncating the corrections. 
The crucial problem for the Bethe approximation today 
is that of finding a small parameter which controls its 
corrections. 

In 1990 Efetov [§} proposed a new formalism for cal- 
culating corrections to the Bethe approximation. Much 
more recently Parisi and Slanina [lfj worked out certain 
details of applying the new scheme to a D-dimensional 
geometry. The advantage of Efetov's formalism is that 
it uses the standard mathematical machinery of saddle 
point approximations, while the geometrical details of 
the Bethe approximation and its corrections have some- 
how been hidden in a clever choice of Lagrangian. These 
authors claimed that the new Efetov Lagrangian: 

• reproduces the Bethe approximation as the zeroth 
order of a saddle point approximation. 

• is able to reproduce the partition function of any 
statistical model with two-body interactions, via 
suitable choices of the Efetov Lagrangian's param- 
eters. For example, the Efetov Lagrangian can 
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reproduce any Potts model, including the Ising 
model, as well as any lattice gauge theory. It is 
important to keep firmly in mind the distinction be- 
tween the Efetov Lagrangian itself and the "target 
theory" - the model which the Efetov Lagrangian is 
being tuned to reproduce, accompanied by a com- 
plete specification of the geometry. 

The second assertion, giving the Efetov theory's physical 
meaning, was based on a saddle point expansion around 
the trivial saddle point. The first assertion, giving the 
connection with the Bethe approximation, was based on 
a saddle point expansion around a different, non-trivial, 
saddle point [251 ]. 

In this paper, section [TT] questions Efetov, Parisi, and 
Slanina's second assertion, that the Efetov Lagrangian 
can replicate the partition function of any target theory. 
It reports that while they are correct that one of the cor- 
rections to the trivial saddle point is the target theory's 
partition function, this correction is only one of many. 
The other (newly discovered) corrections dominate over 
the target partition function, contain two or more spins 
for each spin in the original target theory, and correspond 
to a sum over geometries reminiscent of dynamically tri- 
angulated quantum gravity [TJ E 03 > but without the 
restriction to simplexes having triangular faces. Even 
though the target theory which the Efetov Lagrangian is 
trying to replicate lives on a cubic lattice with D dimen- 
sions and N nodes, the dominant corrections to the trivial 
saddle point live on random geometries with many nodes 
and no specific dimension. The dominant random geome- 
tries are locally tree-like, which perhaps explains why the 
Bethe approximation emerges from the non-trivial saddle 
points. 

This is a very illuminating result. It reveals that im- 
proving on mean field theory in D dimensions is a matter 
of defining an appropriate sum over geometries (prefer- 
ably using an integral formulation) , and analyzing it with 
techniques from statistical mechanics and quantum grav- 
ity. If the relative importance of D-dimensional lattice 
geometries vs. others can be tuned smoothly, then a 
mean field theory with systematic corrections can be de- 
veloped [26| . Assuming that smooth reweighting can ar- 
rive at a regime where tree-like geometries dominate, the 
mean field theory would likely be governed by the Bethe 
approximation. Even if smooth reweighting is not pos- 
sible, any non-trivial saddle points will offer an approx- 
imate mean field theory with systematic corrections as 
long as lattice geometries are dominant. 

As it stands Efetov's theory does not fulfill its promise: 
lattice geometries do not dominate. Moreover, the sum 
over geometries grows factorially with the volume, not 
exponentially, and is therefore in need of regularization. 
This is a familiar problem from quantum gravity and 
string theory [n], [12] , where almost universally the sums 
over geometries with different topologies diverge, and is 
a sign of a nonextensive entropy [161 ]. See [171 ] for a 
recent discussion of the link between nonextensive en- 
tropy, string theory, and replica symmetry breaking in 



frustrated systems. 

Sometimes a failed (divergent) theory is a success, if 
it points the way to other better theories, or gives phys- 
ical insight into a problem - this is the motivation for 
continued research in quantum gravity and strings. The 
Efetov theory is a similar sort of success: it suggests that 
a similar sum over geometries, suitably regularized using 
expertise from quantum gravity and string theory, could 
provide a systematic scheme for improving on the Bethe 
approximation. 

The key technical question is whether a suitable reg- 
ularization can be found. Efetov's theory considers only 
discrete geometries; therefore a lattice cutoff is useless. 
Alternative regularizations like those found in dynami- 
cally triangulated quantum gravity may be useful. Topo- 
logical regularizations [TJ [lj] have found some success 
especially in 1 + 1 dimensions - these restrict the sum over 
geometries to planar geometries (topological spheres with 
no holes), or tori (a single hole), etc. More recently, an- 
other restriction on the sum over geometries - Causal Dy- 
namical Triangulations - has been shown to successfully 
regularize quantum gravity [131 ]. There is good reason 
to hope that physicists in quantum gravity and string 
theory, armed with their expertise in designing and an- 
alyzing sums over geometries, could find a sum over ge- 
ometries suitable for calculating corrections to the Bethe 
approximation. They now have the possibility of making 
a direct and immediate impact on low energy physics, i.e. 
condensed matter theory. 

In other words, the discovery that the Efetov theory is 
a sum over geometries suggests that the Bethe approxi- 
mation - the most powerful mean field theory in existence 
- is best understood from a perspective where geometry 
is not pre-determined or static, but a physical variable. 
From this perspective, the Bethe approximation is the 
low energy limit of some - as yet unknown - sum over 
geometries similar to those found in quantum gravity. 
Perhaps the historical development of the Bethe approx- 
imation via heuristic arguments was both fortuitous and 
accidental; perhaps future physicists will find that rigor- 
ous derivation of the Bethe approximation is an exercise 
in quantum gravity. In my opinion, this new perspective 
is the most important contribution of the present pa- 
per. For a parallel perspective, see Benedetti and Loll's 
recent argument that averaging over geometries may im- 
prove the analysis of critical behavior in D dimensions by 
removing lattice artifacts [14|, [l5j . 

Section IIIII turns to the non-trivial saddle points - the 
ones producing the Bethe approximation and its per- 
turbative corrections. Where Efetov, Parisi, and Slan- 
ina concentrated on developing the perturbation theory 
around the saddle points, here the focus is on a full char- 
acterization of the saddle points themselves. Section IIIII 
begins by showing for the first time how to integrate out 
the exact symmetries of the Efetov Lagrangian, a chal- 
lenge which Parisi and Slanina had noted but not solved. 
It next works out the details of the saddle points and their 
Hessian, giving special attention to a novel mathematical 
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challenge which was not discussed by Efetov, Parisi, and 
Slanina: their theory's Hessian and bare Green's func- 
tion are nonnormal, meaning that [H, H^] ^ 18]. Such 
matrices have distinct sets of right and left eigenvectors, 
and do not obey many of the usual theorems from linear 
algebra. In particular, the right eigenvectors are not or- 
thogonal to each other, and the same applies to the left 
eigenvectors. Nonnormal matrices have been of partic- 
ular scientific interest in describing non-equilibrium and 
chaotic phenomena [HI and also dissipation [2(| , but 
I am not aware of any previous discussion of nonnormal 
bosonic field theory. Lastly, section IIVI summarizes the 
principal contributions of this paper and discusses possi- 
bilities for further research. 



II. THE TRIVIAL SADDLE POINT IS A SUM 
OVER GEOMETRIES 

Efetov, Parisi, and Slanina discuss a class of target 
theories T: 

• The geometry is described by a graph G, with V 
nodes and L links. 

• Degrees of freedom reside on nodes of G. 

• The interactions are restricted to pairwise (two 
body) interactions. 

Efetov, Parisi, and Slanina propose a field theory which 
they claim reproduces the partition function Z T of any 
such target theory T. More specifically, they claim that 
Z T is equal to 

Q = A~ y [ln J d^<hPe Co+XCl - In / d<0#e £o ]. (1) 

If the target theory inhabits a D-dimensional cubic 
lattice, the Efetov Lagrangian [27j is 

£ = £o + \£i= - E i'vds 1 ^vds 2 {l' 1 ) SlS2 

vds\S2 

+ x zZ^ U d^ds^ v _i ds - (2) 

vs 

The complex variables ip are link variables which live at 
each edge on the graph. They are completely specified 
by three indices: the position v, the direction d of the 
link, and the spin index s. These -0 variables are the 
actual degrees of freedom in the Efetov model, and are 
completely distinct from the degrees of freedom of the 
target theory, which are encoded in the sums over the 
spin index s. 7 is a matrix describing the interaction 
between two spin degrees of freedom in the target theory; 
for the Ising model 7 SlS2 = exp(/3siS2). Similarly, £ s 
is a factor encoding the magnetic field. The notation 
v — d used in equation [5] specifies a particular neighbor 
of vertex v. namely the one which is reached by going 
along the d axis in the negative direction. The notation 




FIG. 1: Part a is a simple example target geometry; node 
1 is a box, node 2 a circle. Part b: the two vertices in the 
perturbation theory for this target geometry. 

adopted here assumes that the target theory lives on a 
cubic lattice, but suitable notations can be easily found 
for other target theories. 

In order to understand the physical content of this La- 
grangian, Efetov, Parisi, and Slanina performed a per- 
turbative expansion around the ip = ip = saddle point 
[28| . Parisi and Slanina stated that the perturba- 
tive expansion around the trivial saddle point contains 
a Feynman diagram exactly equal to the partition func- 
tion of the original theory plus diagrams equal to the 
"partition function of the same model wrapped several 
times around." They expect that "in the thermodynamic 
limit the wrapping is unessential" because the wrapped 
geometries differ from the target geometry only at their 
boundaries. 

The main point of this section is to show that the per- 
turbative expansion around the trivial saddle point is ac- 
tually a sum over many geometries differing from the 
target theory throughout their volumes, and to analyze 
the physics of this sum. In fact Q is equal to 

E7 VH 2r (3) 

The first contribution to this sum, , is the desired par- 
tition function of the target theory T. However, there are 
also infinitely many other contributions which cannot be 
neglected: these terms Zj are partition functions of the- 
ories with r times as many vertices as the target theory. 
For instance, if the target theory inhabits a graph with 
100 vertices, then there will be terms with 200, 300, 400, 
etc. vertices. The only way to remove the higher order 
terms is to remove all loops from the graph inhabited by 
the target theory, thus restricting oneself to working on 
trees. 

To show that Q contains a sum over geometries, let 
us work out the details of a particularly simple target 
theory living on two nodes connected by three links, as 
shown in figure [T^,. In this case the Efetov Lagrangian is 

3 

£ = £ + A£i= - EE^WT- 1 )^ (4) 

d=l s±s 2 
s 
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FIG. 2: The second order geometries for our example. Part 
a - the spanning trees. Part b - the disconnected geometry, 
weight 1/2. Part c - the connected geometry, weight 3/2. 
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FIG. 3: The third order geometries for our example. The 
weights are 3 (a), 1 (b), and 1/3 (c). 



In the perturbative expansion C\ generates two vertices, 
each with three legs. As illustrated in figure [TJd, the legs 
are labeled; the first vertex has legs labeled a, b, and c 
while the second vertex has legs labeled A, B, and C. The 
Efetov Lagrangian's bare propagator connects a with A, 
b with B, or c with C. In the formulas corresponding to 
each Feynman diagram, one puts a factor of 7 SlS2 for each 
instance of the propagator, a A ^ s for each instance 
of the first vertex, and a A ^ s £| for each instance of the 
second vertex. 

Zj , the first order contribution to equation [TJ con- 
tains one instance of each vertex, and there is only one 
way to connect the two, shown in figure [TJa,. This graph 
is exactly the same as the target theory's geometry, so 
this first order correction produces exactly the partition 
function of the target theory. However the second order 
correction contains two instances of each vertex, and 
now one must think a bit. There are (2!) possible ways 
of connecting them; one factor of 2! can immediately be 
removed by constructing spanning trees, as shown in fig- 
ure [2ji. We now have (2!) 2 possible choices about how 
to connect the four remaining pairs of legs; figures ^Bp 
and[2t show the two second order diagrams. The first is 
a disconnected diagram with symmetry factor 1/2, and 
is removed by the logarithm in equation [1] The second 
diagram connects the two spanning trees, contains four 
nodes instead of the two nodes present in the original 
target theory, and has symmetry factor 3/2. 

The third order contribution Zj can be found in a 
similar fashion, by first constructing the spanning trees, 
and then enumerating the (3!) 2 possible ways of connect- 



ing the six remaining pairs of legs. There are of course 
disconnected geometries that are again removed by the 
logarithm. However there are also three topologically 
distinct fully connected geometries, which are shown in 
figure [3] along with their symmetry factors. The best 
way to figure out the geometries and symmetry factors 
reliably is to do the enumeration by hand. 

Consider now evaluating Zj for a theory inhabiting a 
more complicated target graph with V nodes and L links. 
Again we will start with r V vertices. Immediately we 
will connect the vertices into r separate spanning trees. 
Each spanning tree will contain V — 1 links, and will 
contain also 2(L — V + 1) dangling legs which still need 
to be connected. However, if the original target theory 
inhabits a tree, then 2(L — V + 1) = and the pertur- 
bation theory stops here; all higher order diagrams are 
disconnected spanning trees and are deleted by the log- 
arithm in equation [TJ In this case Q is exactly equal to 
the partition function of the target theory, precisely as 
predicted by Efetov, Parisi, and Slanina. 

However, if the target theory does not inhabit a tree, 
there will be an additional phase of connecting the re- 
maining dangling legs, and thus tying the spanning trees 
together. This results in many other Feynman diagrams, 
each corresponding to an alternative geometry, a different 
way that the r V vertices could be connected together. 
The logarithm in equation [TJ removes disconnected ge- 
ometries. We now turn to an analysis of the resulting 
sum. 



A. The dominant geometries 

Which geometries dominate the sum? This question is 
answered by summing their symmetry factors. Section 
III A II argues that the sum of the symmetry factors of all 
fully connected geometries of order r is approximately 
(rl) ~ V . If the target theory inhabits a D-dimensional 
lattice then L — V = (d—l)V, so this sum is rather 
large. For instance, on a 10x10 lattice the second con- 
tribution to Q, Z%, will be roughly (2!) 100 times as big 
as Zi , the partition function of the target theory. This 
large factor is somewhat counterbalanced if the model 
has a positive free energy density imposing a cost for 
adding nodes. Nonetheless, because there are geometries 
with r arbitrarily large, eventually the (r\) L V growth 
will overwhelm the free energy density. Put another way, 
the nonextensive growth of the entropy overrides the ex- 
tensive growth of the free energy. In summary, large r 
geometries completely dominate the sum. 

The dominant geometries are locally tree-like, in the 
sense that the only non-backtracking walks which return 
to their starting point are infinitely long. This is obvi- 
ous from section [ITjs discussion of how to enumerate ge- 
ometries, by first constructing r spanning trees and then 
connecting them in all possible ways. Each node lies on 
one of the r different spanning trees. In order for a non- 
backtracking walk to come back to its origin, it must 
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necessarily step off of the origin's spanning tree. Once 
the walk is on a different tree, each subsequent step will 
have a probability of order r _1 of returning to the ori- 
gin's spanning tree. Since large r geometries dominate, 
the average loop length is infinitely large. 

The Bethe approximation is exact on trees, and a good 
first approximation on geometries which are locally tree- 
like. Perhaps this explains, on a physical level, why the 
Efetov Lagrangian exhibits saddle points related to the 
Bethe approximation. 



1. Estimate of the connected symmetry factors 

First of all, the nodes in a geometry are generated by 
an exponential e XCl ; the r instances of each node in the 
target geometry are accompanied by a factor of (?*!) . 
Connecting up all the nodes requires rL links, giving a 
factor of (r\) L . Therefore the sum of all the order r sym- 
metry factors (of both connected and disconnected ge- 
ometries) is (r!) 

I now argue that the proportion of disconnected ge- 
ometries in this sum is very small. I write C r for the 
sum of the symmetry factors of the order r connected 
geometries, and D r for the sum of the disconnected sym- 
metry factors. Thus C r = {r\) L V — D r . Now the 
number of disconnected geometries D r is simply related 
to the numbers of connected geometries at lower or- 
ders: D l = 0, D 2 = \{C x f , D 3 = dC 2 + i(d) 3 , 

D 4 = CiC 3 + i(Ci) 2 C 2 + i(Ci) 4 + i(C 2 ) 2 , et cetera. 
I am not aware of any way to write a closed form expres- 
sion for the total value of D r except in the special case 
L = V, but its largest single part is the C\C r -\ term, 

which is smaller than (r!) L V by a factor of (i) 

All other contributions are much much smaller, of order 

( r ( r _i) ) or l ess - Therefore I conclude that if L — V is 
large then the connected weights sum to approximately 

i \\ L ~ V 
(r!) 



2. Divergence of Efetov' s theory 

The perturbative expansion around the trivial saddle 
point diverges. I have already argued that if L — V is 
large then the total number of order-r geometries is ap- 
proximately (H) L . Because the free energy density of 
each geometry is bounded from above by a geometry- 
independent value /, the partition function of each ge- 
ometry is bounded from below by exp (— (3rVf). Since 
the number of geometries (r\) L V increases faster than 
exponentially and each geometry gives a positive contri- 
bution, the terms Zj w X rV (r\) L ~ V exp (— flrVf) in the 
perturbative expansion around the trivial saddle point 
increase without bound, and the entire series diverges. 
The fact that the individual terms in the series are all 
positive provides strong evidence that not only the per- 



turbative expansion of Efetov's theory diverges but also 
the theory itself. 

This divergence has nothing to do with the volume: it 
occurs even when the target geometry has only a small 
number L of links, in which case Efetov's path integral 
integrates over only a finite number of degrees of free- 
dom. Instead the divergence is a result of interaction 
terms in the Efetov lagrangian, much like the way that 
(j) 4 theory diverges if the interaction has the wrong sign. 
Unlike 4 theory, the interactions in Efetov's theory are 
complicated enough that the divergence is not easily vis- 
ible unless one expands perturbatively in the interaction 
constant and obtains a sum over geometries. As we have 
seen, the perturbative picture of the divergence is that 
the number of geometries grows factorially, while their 
free energies grow only exponentially. 

Within the framework of the Efetov Lagrangian, the 
only way to control the divergence is by treating target 
geometries which have less than two loops or by setting 
A = 0. Quantum gravity manifests a similar divergence, 
but it has been regularized by constrainin g th e sum over 
geometries, either by fixing the topology [ill [l2| or via 
Causal Dynamical Triangulations [13| . 

B. Is the geometrical interpretation valid? 

We have just established that the perturbative expan- 
sion around the trivial saddle point diverges, which is 
standard, since as a rule the perturbative expansions 
of non-trivial field theories are only asymptotic. To be 
more precise about asymptotic theories, their successive 
perturbative corrections first decrease term by term but 
eventually start increasing factorially, due to the pertur- 
bative expansion's neglect of non-analytic features of the 
original path integral. Nonetheless it is well known that 
the first few terms of an asymptotic expansion can give 
a very good approximation to the correct result. 

There are several reasons to think that Efetov the- 
ory's divergence is much more malignant, and to doubt 
altogether the validity of perturbation theory around the 
trivial saddle point: 

• Unlike an asymptotic theory where the coefficients 
of the perturbative contributions decrease for a 
while before starting to diverge, this series ex- 
hibits exponentially large coefficients already at 
the second order. The second order contribu- 
tion to the trivial saddle point is approximately 
(A exp (— Pf)) V (2\) L ~ V bigger than the first order 
contribution, so that the small parameter A is coun- 
terbalanced by a large geometrical factor. The free 
energy term exp (—/?/) is freely adjustable by addi- 
tive contributions to the Lagrangian, and therefore 
should not be regarded as distinct from the pertur- 
bative parameter A. 

• The coefficients increase not only factorially, but as 
a factorial taken to a very large power: (r\) L V . 
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• Every order of the expansion contributes with the 
same sign. 

The geometrical interpretation of the Efetov Lagrangian 
- including both its link with the target theory T and its 
interpretation as a sum over geometries - rests on the va- 
lidity of this perturbative expansion. Since this is ques- 
tionable, the geometrical interpretation must be taken 
with a grain of salt. Probably the only way to rigorously 
justify a connection between the Efetov Lagrangian and 
any target theory will be to introduce some type of reg- 
ularization to the Efetov Lagrangian. 

III. THE NON-TRIVIAL SADDLE POINTS 

The previous section discussed the physical meaning of 
Efetov's theory which is based on the trivial saddle point. 
This section turns to the non-trivial saddle points, which 
are related to the Bethe approximation. Efetov, Parisi, 
and Slanina have already devoted much ink to develop- 
ing the perturbation theory, in hopes that it will provide 
systematic corrections to the Bethe approximation. The 
focus here is on giving a full account of the saddle point 
structure, including Goldstone bosons, complete enumer- 
ation of the saddle points, computation of the Hessian 
and its eigenvalues, and the Hessian's nonnormality. 

A. U(l) symmetries and Goldstone bosons 

The Efetov Lagrangian, given in equation^ possesses 
many exact U(l) symmetries. Each of the links can be 
multiplied by an overall phase: VWs - ► ^vds^ x ' d ■ There 
are L links, but only V different linear combinations of 
the L phases actually appear in the Efetov Lagrangian, 
one for each term in C\\ these linear combinations are the 
sums 9 V = J^d'fivd ~ ^v-dd ^ich obey the constraint 
J2 V 9v = 0. Therefore the Efetov Lagrangian has L — 
V + 1 exact symmetries. There is a strict analogy to the 
following integral with a single exact £7(1) symmetry: 

J ##e- m ^ +A ^ 2 . (5) 

The exact symmetries create difficulties for the evalua- 
tion of non-trivial saddle points. Goldstone bosons arise 
from each of them due to spontaneous symmetry break- 
ing. Goldstone bosons, especially ones caused by exact 
symmetries, are invitations to integrate exactly, sepa- 
rately from any saddle point approximation. Indeed, any 
attempt to include a Goldstone boson in a saddle point 
approximation will give an infinite result. 

B. Exact integration 

The exact U(l) symmetries may be integrated as fol- 
lows: 



• Factor out a phase e 1 ^ from each link. One does this 
by changing to angular variables ip v d s — r v dsS l ^ vdB , 
adding <fi v d,s=o to the other phases {4> v d, s^o — * 
<ftvd,s=£o + <Pvd,s=o), an d factoring out the common 
phase <j> v d,s=o- A Jacobian term Y^vds m IVw«l 
must be added to the Lagrangian. 

• Notice that the Lagrangian depends only on the lin- 
ear combinations 9 V = X^0„<i, a =o - <t> v -d,d,s=w 
and therefore integrate out the L — V + 1 other 
linear combinations, which just multiplies Efetov's 
path integral by a constant Jacobian J — 1 and by 
(2vr) L - y+1 . 

• Impose the constraint ®v = by inserting a 
delta function into Efetov's path integral. Without 
this constraint the 9 integrations would decouple 
and the path integral would be exactly zero. 

The final Lagrangian is 

c= - ^'i S1 tp vdS2 ('y' 1 ) slS2 +^2^\^vd S \ 

vds\S2 vds 

+ \Y,te" v Wvdrfv-t*,- (6) 

vs 

The tjj's are in radial coordinates ipvds = r v d s e l ^ , ' d3 , and 
4>vd,s=o = 0. There is some liberty to convert the i^vd,s=jto 
components back to Cartesian coordinates. This sim- 
plifies the perturbation theory but mutilates the Efetov 
theory's symmetry with respect to the spin index [29| . 

The presence in the action of e lBv , a complex tran- 
scendental function, raises some concern about apply- 
ing the saddle point approximation to the 9 integra- 
tions. Unfortunately exact integration of the 9 vari- 
ables would add to the Lagrangian a transcendental term 
of the form J2 V ln /(^ S s 6 n <zV>i>dsVVd,ds)' prohibiting 
further progress. 

C. The saddle points 

As is usual when doing mean field theory, we look for 
a spatially uniform ground state; one that depends on 
neither the site index v nor the direction index d. It 
is helpful to first consider the saddle point equation ob- 
tained from the final Lagrangian given in equation [5] but 
ignoring the logarithm. We take the first derivative of 
the Lagrangian and obtain 

The saddle point solutions scale with A as \^ 1 /( 2D ^ 2 '> ■ if 
we remove A and constrain ip to be real then the saddle 
point equation becomes the Bethe-Peierls approximation: 

(8) 

Sl 
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I have set 9 to zero because this is the only spatially 
uniform value which will satisfy the global constraint. 
On the other hand, the Lagrangian's first derivative with 
respect to 9 V is «A^ s £ s V>s D ; not only non-zero but also 
imaginary. The real part of the first derivative vanishes 
at the non-trivial saddle points, but the imaginary part 
does not. Perhaps this is OK, since the imaginary part 
contributes only a phase. As a consequence of this diffi- 
culty, 9 will have a vacuum expectation value. However, 
the global constraint on 9 removes the zero-momentum 9 
mode from the Hessian, so the vacuum expectation value 
has no influence on the perturbation theory. 

Equation [7] is almost the same as the saddle point 
equation which is obtained from the original Efetov La- 
grangian in equation [5] The difference is that the sad- 
dle point equation before U(l) integration has a U(l) 
symmetry which must be spontaneously broken by hand, 
while after integration the phase of ip s =o is constrained 
to be zero and the U(l) symmetry is no longer present. 
Other than this procedural detail, the two equations, and 
their solutions, are exactly the same. 

Equation [7] contains 25* equalities which must be sat- 
isfied, where S is the number of values that the spin can 
take. Since the imaginary part of ip s= o is zero, there are 
only 2S — 1 unknowns; one should not expect complex 
solutions. On the other hand, real solutions are possible, 
since in this case there arc only S equalities and S un- 
knowns. In particular, each of ip's components, with the 
exception of ^> s =o, may have either a negative or positive 
value. Some solutions of the Bethe approximation actu- 
ally do have negative values, even though this conflicts 
with the physical reasoning leading up to the Bethe ap- 
proximation, in which ip a is understood as the value of a 
partition function. 

When the target theory is the Ising model, the Bethe 
approximation has more than one solution. I discuss 
only the ferromagnetic regime where [3 J > 0. Two non- 
zero real solutions are present at all temperatures, one 
of which gives a negative value to ip s= i- Below a critical 
temperature two more solutions appear. As a general re- 
sult for all target theories on cubic lattices, these Bethe 
saddle points are not local maxima of the Lagrangian. 
This may be verified by examining the Lagrangian's Hes- 
sian (see section IIII D[) at the Bethe saddle points and 
noting that it has a positive eigenvalue corresponding 
to multiplications of the saddle point solution ip by a 
real constant. As a consequence, the states correspond- 
ing to the Bethe saddle points are not stable; they are 
not true ground states, and it is impossible to define a 
valid perturbation theory around them. Efetov's original 
program for calculating corrections to the Bethe approx- 
imation depended crucially on doing perturbation theory 
around the Bethe saddle points; the non-existence of such 
a theory means that the Efetov Lagrangian is useless for 
Efetov's original goal, though perhaps quite useful for 
showing a way forward via connections with quantum 
gravity. 



1. Additional saddle points resulting from the logarithm 

If one does not ignore the logarithm in the final La- 
grangian shown in equation [5J the resulting saddle point 
equation is not quite the Bethe-Peierls approximation 
any more: 

= £ ls 2Sl (xu ft?- 1 + \^). (9) 

si 

When X 1 /^ ^ 1 ) is very small, equation [9] exhibits all 
the Bethe solutions from equation [SJ plus four more in 
the particular case of the Ising model. As \ 1 /( D ~ 1 '> is 
increased, the saddle points change continuously and no 
new ones appear. The stability of each saddle point re- 
mains the same as well, except at large A. Finally at 
A ~ 0.1 all of the solutions disappear. The Bethe saddle 
points all scale as X^ 1 /^ 20 ^ 2 ) ; so at these saddle points 
the last term in equation [5] is of order X 1 /^ ^ 1 ) and can 
be neglected, and the corresponding term in the Hessian 
is also very small. Therefore all the Bethe saddle points 
are unstable. 

We turn to the other four saddle points, in the limit 
of small X 1 ^ ^ 1 ^. The extra saddle points all disrupt 
the balance between the two terms in the original Bethe 
approximation, so that either ij) s or X£ s i/; 20-1 is negli- 
gible. As a result these saddle points have no relation 
to the Bethe approximation. Two of them have ip s =o = 
±&=i = v /(e /3J ±e- /3J )/2- These zero- magnetization 
solutions are notable both because they don't depend on 
the magnetic field at all and because the + solution is 
stable. The — solution is not stable. There are also two 
ferromagnetic saddle points where one component scales 
as X- 1 ^ 20 - 2 ^ and the other scales as A 1/(2D ~ 2) . These 
are both unstable. 

One intriguing aspect of saddle points is their depen- 
dence on coordinate transformations. The extra term 
in equation [5] is entirely due to a transformation from 
Cartesian coordinates to angular coordinates. The same 
term could be removed again by transforming from r to 
w = r 2 . If one wanted to preserve the Bethe approx- 
imation while playing such games, one would need to 
introduce a term which would not disturb the A -1 ^- 0-1 ) 
scaling of the solutions. Therefore the coordinate trans- 
formation would have to depend on A. 

D. The hessian 

Evaluation of a saddle point and its corrections starts 
with computing the Hessian (second derivative) of the 
Lagrangian at that saddle point. The inverse of the Hes- 
sian will be the bare Green's function used in pertur- 
bative corrections to the saddle point, while a Gaussian 
integral incorporating the Hessian will determine the sad- 
dle point's weight. 

Because we have removed the £7(1) symmetries, we 
begin using the following degrees of freedom explic- 
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itly: if) = (r,4>,6). The following rescaling simplifies 
the Lagrangian's form around the saddle point: ip s = 
+ r s /p s ) exp (t(p s /p s ), where p s = £« We also 

rescale 8 — ► 9/\p\. In these rescaled variables, the La- 
grangian is 

- ^i^2(7 _1 ) SlS2 (l + r v dsjp si ) (1 +r vdS2 /p S2 ) 

vdsiS2 

X exp (l(j) vdsi /Psi - rfvdsa I Ps 2 ) 

+ XJ2p 2 s e l6Jm n d (1 + r vds / Ps ) (1 + r v _ id Jp s ) 



X exp (t(j) v ds/ps ~ l^y-ldsl Ps) 

+ ^ (ln^s +ln|l + r vds /p s \). 



(10) 



The Hessian is found by extracting all the second order 
terms from the Lagrangian, and is 

jjrr _ _ 2A5 VlVa 5d 1 d 2 + XS SlS2 S VlV2 (l — Sd 1 d 2 ) 

~~ Psx &s 1 s 2 &v 1 v 2 5d 1 d 2 + X5 SlS2 (S T + S + 6*2), 

H U = -2PA(l-8 slS2 )P5 VlV J dld2 

— 2(\ — a sl )PS slS2 S vlV2 5d 1 d 2 

- XPS mV2 (l-S dld2 ) + XP(S T + S-S 2 ), 

Tree _ _ X c 

n ^0 VlV2 j 

jjr<t> _ (TJ<l>r\ T 

= i\P6 VlV2 (l - S dld2 ) + *\P(S T -S- S 2 ), 



H 



Or 



(H rb ) =i\p s (5 VlV2 +S). 



(11) 



In the above equations <5 Sl s 2 is the Kro- 
necker delta function, S is the stepping matrix 



Svi di , 



v 2 d 2 



= 5 



trix S2,v 1 d\,v 2 d 2 



v\,v 2 +d 2 



, ^2 is the two-step ma- 
(1 - ^1^2)^+^, 



a 1/ % 2 1/2 ^; D ^ 2 - D (i- 1 ) slS2 , 



v 2 +d 2 ' 



A is 



a is 



A - 

■"-S1S2 — 

o.x = ^Ir^E^V^^ 1 )^, and P is the 
projection operator P SlS2 = <5 SlS2 (l — <5osi)- One must 
keep in mind the global constraint on 6, and also that 
A, a, and p all scale with various powers of A. 

In the zero-momentum sector the Hessian is consider- 
ably simplified: 

H rr = -2A6 dld2 +2X6 slS2 (2-6 dld2 ) 

~ Psi 3s 1 s 2 3did2> 

H H = - 2 PA(l-5 SlS2 )P8 dld2 

+ 2a Sl P5 SlS2 5d 1 d 2 - (12) 

E. The hessian is nonnormal and pseudo-hermitian 

The degrees of freedom break into two classifications: 
<j> and 0, which are always accompanied by a factor of 1, 



and r which occurs alone. Broken into these sectors, the 
Hessian has the form 



H = 



X 



a x 
T (3 



(13) 



where a and (3 are real and hcrmitian and \ is imaginary 
but not hermitian. Therefore, H 's determinant is always 
real, and H itself is neither hermitian nor anti-hermitian. 

Any matrix in this form is also nonnormal unless ax — 
X/3, which translates to H^H^ + H r6 H e ^ = H^H^ 
for the case at hand [3C} • This equality is not verified, so 
in fact the Hessian is nonnormal at the non-trivial saddle 
points [311 ]. 

The Efetov Lagrangian's nonnormality results from 
treating the underlying geometry as a directed rather 
than undirected graph. Any theory without a simple 
symmetry between tp and ip will have a nonnormal Hes- 
sian. Even a free field theory C = —tpDip is nonnormal 
if [D, D T ] ^ 0. 

Within the class of nonnormal matrices, the Hessian is 
specifically a pseudo-hermitian matrix [321 ] , which means 
that the hermitian conjugate of H has a special rela- 
tion to H itself: H = CH^C^ 1 . In our case C is the 
1 



conjugation matrix 



-1 



As a consequence of the 



pseudo-hermitian property, if the set of eigenvectors of H 
is complete then all the eigenvalues of H either are real or 
occur in complex conjugate pairs, and C maps the right 
eigenvectors to the left eigenvectors j20| . 

If the target theory possesses translational invariance, 
then the Efetov Lagrangian shares the same symmetry. 
In this case the Hessian commutes with the translation 
matrix, and decouples into sectors corresponding to each 
value of the momentum. Because the Hessian's nonnor- 
mality is a symptom of spatial structure, the Hessian's 
zero- momentum sector is normal, as seen in equation ll2l 



1. The nonnormal saddle point integral 

When using the saddle point approximation, the 
weight of each saddle point is determined by a product of 
several quantities, among which is the Gaussian integral 



dip exp (—tpHif). 



(14) 



Switching to a momentum basis, one obtains a product 
of Gaussian integrals which have exactly the same form 
as integral 1141 The zero-momentum integral is normal, 
but the other integrals are nonnormal. 

I am not aware of any published way to do this in- 
tegral exactly when, as in our case, H is nonnormal. 
The standard results for Gaussian integrals cover only 
cases where either ip is a Grassman variable or else ip 
is a scalar variable and H is a normal matrix. In the 
second case a suitable change of variables will factor- 
ize integral [TJ] into a multiple of many one-dimensional 



9 



Gaussian integrals. The integral is then equal to 

0) S/2 exp (-±Tr (In -#)) = (27r) B/2 (det {-H)y 1/2 
if H is negative definite, and diverges otherwise. In con- 
trast, because our Hessian is nonnormal, its eigenvectors 
are not orthogonal, and there is no change of variables 
that will cause integral [14] to factorize. 

Theorem: If H is symmetric and negative definite and 
its eigenvectors are complete and contained in the matrix 
W, then 



dip exp(-^Htp) = (-l)V' 



(2tt 



B/2 



det(~H 



1/2 ■ 



(15) 



where Wi are the eigenvalues of W, a = 
\ J2i (sign(wi) - l) 2 , and uj = £\ phase(wi) . 

Proof: H has the decomposition H = WXW -1 , 
where A is the diagonal matrix composed of if's eigen- 
values. We change coordinates ip ~~ * Wip, obtaining 
(-lfdet [W)J dip exp^ifw^WXip). Now imagine do- 
ing the individual integrals one by one, completing the 
square each time. This process of completing the squares 
is mathematically equivalent to computing the LU de- 
composition of W^VKA. LU decomposition means fac- 
toring a matrix into two matrices L and U, where L 
is zero above the diagonal, U is zero below the diag- 
onal, and the diagonal elements of L are equal to one 
[Tl |. If the LU decomposition exists then it is unique. 
After completing the squares, the z-th integration con- 
verges if Re(Uu) < 0. The result of all the integra- 
tions is (-l) CT det (W) (2ir) B/2 (TL, -Uu)~ 1/2 . The prod- 
uct of — E/'s diagonal elements is just the determinant of 
— W^VKA, giving formula ITSI 

The last step is to establish that the LU decomposi- 
tion exists and that Re(Ua) < by constructing L and U 
more explicitly. W^W is positive definite and therefore 
has a Cholesky decomposition W^W = G^G, where G is 
zero below the diagonal and its diagonal elements are real 
and positive [lfj. Defining the matrix gij = SijG^ 1 , we 
obtain W^WX = G^gg^GX. Therefore the LU decom- 
position exists: L = G^g and U = g~ 1 GX. The diagonal 
elements of U are G 2 { Xi, so the integral exists if H is 
negative definite. Q.E.D. 

Lemma: If H has the form given in equation 1131 is 
negative definite, and its eigenvectors are complete, then 



(16) 



In this case the Gaussian integral must be real because 
J dtp is symmetric under the transformation ip — > Cip, 
which reverses the sign of the imaginary part of ipHip. 
The only possible values for the phase are ±1. If H is 
normal then the ambiguity in sign is correctly resolved 

1/2 

by taking the absolute value of det(H) 1 , and it seems 
likely that the same applies to nonnormal H. 

This theorem does not address the possibility that 
equation [T?] might converge even if one or more eigen- 
values of H is zero or positive, since it does not exclude 



the possibility of using a set of coordinates other than 
tp — > Wtp. If H is in the form of equation [T31 then con- 
vergence is assured as long as a and (3 are both negative 
definite. Whether violation of this condition ensures di- 
vergence remains to be seen. Of course divergence is 
assured if the zero-momentum Hessian is not negative 
definite. 

One can pursue a perturbative strategy, dividing iJ's 
real part R which is hermitian from its imaginary part /, 
and treating / as a perturbation: 



H = 



X. 



a \ 
T (3 



= R + I, 



R 



a 




X 



o X 
T 



(17) 



dip exp (-ipHip) 



exp (7--^,- / • -i) f dip exp (^-ip Rip + ip ■ I) 
2 dl dl J 2 



(18) 



1=0 



Clearly the remaining integral diverges unless R is nega- 
tive definite. Summing all orders in perturbation theory, 
we obtain 

(2^) S/2 exp {-\Tr [In -R] ) (19) 



x exp 



» (2tt) B/2 exp ( - -Tr [In -H] ). 

This is exactly the same result that one obtains rigor- 
ously via diagonalization when H is normal and negative 
definite, except that all non-perturbative information has 
been lost. The perturbation theory suggests that conver- 
gence depends on the spectrum of R, while we know that 
if H is normal then convergence depends on H, not R. If, 
as in our case, H is not normal, then perturbation theory 
gives no indication about whether there might be finite 
non-perturbative corrections to equation 1191 or even of 
the conditions for divergence or convergence. There is 
reason to be suspicious of equation 1191 since at most of 
the Efetov theory's saddle points the perturbation is not 
small. 



F. Spectra 

As I have already mentioned, all but one of the Efe- 
tov theory's non-trivial saddle points have at least one 
positive eigenvalue in their zero-momentum sector and 
are therefore guaranteed to diverge. The single negative 
definite saddle point has a magnetization of exactly zero, 
does not respond to an external magnetic field, and is 
not a solution of the Bethe approximation. This saddle 
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point is negative definite at all momenta, ensuring that 
integral Q3] converges. 

The spectra of the saddle points which satisfy the 
Bethe equation all scale linearly with A. Therefore the 
bare Green's function at the Bethe saddle points scales 
as A -1 . Since the perturbative vertices are proportional 
to A, any perturbative diagram with more lines than ver- 
tices will diverge proportionally to a negative power of 
A. 

There are two saddle points which are constant with 
respect to A, including the stable saddle point. Their 
spectra are constant with respect to A, as are the bare 
Green's functions. Therefore perturbative corrections are 
controlled by the small parameter A, with one power of 
A for each vertex. 

There are also saddle points where one component 
scales as X^ 1 /^ 20 ^ 2 ) and the other scales as X 1 /^ 20 ^ 2 ). 
The scaling behavior of their spectra is harder to an- 
alyze, but there are eigenvalues in the zero-momentum 
sector which scale with X~ D /( D ~ 1 \ as well as ones that 
scale with X ^ - 1 ^. 



G. The free energy density 

The free energy density, neglecting both constants and 
perturbative corrections, and assuming that the Hessian 
is negative definite, is 

/3f = ln|pl- J Dln|p |+2£»^ln| /9s | 

S 

+ ^Tr\n{-H) + {D-l)X\p\ 2 . (20) 

If the saddle point is dominated by a specific number 
of vertices, then the partition function is equal to A to 
some power and the free energy density / is proportional 
to log A. This is the case for the unconditionally stable 
saddle point, where p and H are proportional to 1, which 
means that it does not involve any vertices at all, much 
like the trivial saddle point. On the other hand, the free 
energy density at the Bethe saddle points, if the Hessian 
were negative definite, would be dominated by the last 
term in equation [20l This term scales as A -1 ^ 15-1 ), so 
these saddle points are superpositions of many different 
geometries. 



H. Reliability of the computations 

Some of the results reported here were obtained non- 
analytically, through a code which computes the saddle 
points and their Hessians and eigenvalues. The temper- 
ature, magnetic field, coupling constants, number of di- 
mensions, and lattice size can all be specified. Many au- 
tomated tests are built into the code, including compu- 
tation of the Hessian three different ways. If the reader 
wants to reproduce or check these results, the code is 



available under the GNU public license and may be down- 
loaded from my web site, www.sacksteder.com. 

IV. CONCLUSIONS AND FURTHER 
CHALLENGES 

The most exciting result of this paper is that sums over 
geometries may be the appropriate tool for understand- 
ing and improving mean field theory in D dimensions, 
and in particular for developing systematic corrections 
to the Bethe approximation. There is considerable room 
for designing various sums over geometries in which the 
lattice geometries are already dominant or else can be 
smoothly reweighted to be dominant. If possible, the 
sums should be finite, be formulated as integrals, and 
have non-trivial saddle points which are local maxima 
of the Lagrangian. In this case they would define in- 
teresting mean field theories with systematic corrections. 
In order to obtain corrections to the Bethe approxima- 
tion, one would hope to find an ensemble which allows a 
smooth reweighting, without phase transitions, between 
lattice and tree-like geometries. Quantum gravity and 
string theory may provide exactly the expertise needed 
to construct the needed ensemble. 

The Efetov theory is, to my knowledge, the first ex- 
ample of a bosonic theory with a nonnormal Hessian and 
Green's function. Nonnormal field theories hold some 
promise for studying non-ergodic and non-equilibrium 
behavior [l9l |. and for the study of dissipation in quan- 
tum systems [13, l23| . This paper proves for the first 
time the value of Gaussian integrals with nonnormal ker- 
nels, which is a prerequisite for the development of saddle 
point approximations and perturbation theory. It also 
opens a host of questions: 

• What are the necessary and sufficient conditions for 
convergence of a Gaussian integral with a nonnor- 
mal kernel? Can the phases and signs in equations 
[T5l and [16] be simplified? Integrals with kernels in 
the form given by equation [T3] are particularly im- 
portant because this is the most general form for 
bosonic Lagrangians that contain no complex con- 
stants. 

• It is well known that nonnormal matrices are not 
fully characterized by their eigenvalues. For in- 
stance, their spectra can be exceedingly sensitive 
to small perturbations, and their powers may show 
a transient behavior which grows much faster than 
powers of the largest eigenvalue. Pseudospec- 
tra, plots of how the spectrum would vary under 
small perturbations of the matrix, are re cog nized 
as showing much additional information [19]. Do 
these issues manifest themselves in the structure 
of perturbative corrections to a nonnormal saddle 
point? 

• What is the physical meaning of Feynman dia- 
grams in a nonnormal theory where propagators 
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are extremely sensitive to perturbations? Does self- 
energy still have a meaning? 

This paper also contributes a more detailed analysis 
of the Efetov theory, culminating in the discovery that 
it is actually a sum over geometries, and also in strong 
evidence that it is divergent and in need of regular- 
ization, probably similar to regularizations in quantum 
gravity. Even without regularization, its interpretation 
as a sum of tree-like geometries argues that further at- 
tention should be given to the 9 integrations. If some way 
were found to do these integrals non-perturbatively, per- 
haps the Efetov theory would yield saddle points which 
are more satisfactory. 
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equations and the cavity approximation |l| . It also is the 
foundation of the belief propagation and message passing 
algorithms which have proved useful in coding theory and 
in combinatorial optimization problems [2ll.[22|. 

[25] There was an implicit assumption, which will not be dis- 
cussed further here, that the two saddle point expansions 
must converge to the same value. 

[26] For theories with a Lagrangian structure similar to the 
Efetov theory, it may be possible to force lattice domi- 
nance by adding a term to the action proportional to the 
number of short self-avoiding loops. 

[27] Efetov, Parisi, and Slanina do not include the perturba- 
tive parameter A in their theory, which is equivalent to 
setting it equal to 1. 

[28] Note that there is no such saddle point if the target the- 
ory inhabits a graph which contains dangling nodes; how- 
ever it is possible that such theories have a saddle point 
solution with values exponentially close to on nodes 
which are far from the dangling nodes. It is also worth 
noting that the trivial saddle point is not a local max- 
imum of the action unless 7 is positive definite. In the 
Ising model this amounts to requiring a ferromagnetic 
coupling. 

[29] Rather than integrating the U(l) symmetries exactly, 
one could adopt a different approach: simply remove 
the Goldstone bosons from the perturbative corrections 
to the saddle points, by truncating the zero sector of 
the bare Green's function. This is equivalent to insert- 
ing L — V + 1 delta functions into Efetov's path integral, 
given in equation [1] and is exactly analogous to inserting 
a delta function S(Im tp) into equation transforming 
it to / dip e ~^ 2 +^ 4 . Clearly this is an effective way 
of removing exact symmetries, although it destroys any 
rigorous connection with the physical content of Efetov's 
theory, discussed in section [TTJ In this approach, the ex- 
tra ip~ term is not added to the saddle point equation. 
Because Efetov, Parisi, and Slanina do not mention this 
term, it would appear that they truncated the Goldstone 
bosons. However, this is not explicitly mentioned, and 
Parisi and Slanina instead write of a cancellation of di- 
vergences at the one-loop level. Their neglect of the Hes- 
sian's nonnormality, which is discussed in section IIII El 
confuses the question still more. 

[30] The mapping from the Lagrangian to the Hessian is not 
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unique: it is possible to add any antisymmetric matrix to 
a and (3, and any matrix to \ as l° n g as it is compensated 
by a matching subtraction from x T ■ This seems to offer 
hope of finding an alternate Hessian which is normal. 
The system of equations for determining such a Hessian 
turns out be quadratic in the unknown matrices, with 
many linear constraints. If any solution exists, deriving 
it is not likely to be trivial. 
[31] This nonnormality was overlooked in the calculations by 
Efetov, Parisi, and Slanina, and their subsequent results 



about the Green's function and about observables beyond 
the Bethe approximation may be incorrect. In particular, 
the expansion of the bare Green's function in terms of 
its eigenvalues and eigenvectors must make an explicit 
distinction between right and left eigenvectors. 
[32] Pseudo-hermitian matrices have recently been of much 
interest in the study of PT symmetric quantum mechan- 
ics [20I . |23| , which is an interesting way of including dis- 
sipation effects into the standard structure of quantum 
mechanics. 



